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Measurement of the Positive Muon Anomalous Magnetic Moment to 0.46 ppm 
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We present the first results of the Fermilab National Accelerator Laboratory (FNAL) Muon g — 2 
Experiment for the positive muon magnetic anomaly a, = (g, — 2)/2. The anomaly is determined from the 
precision measurements of two angular frequencies. Intensity variation of high-energy positrons from 
muon decays directly encodes the difference frequency w, between the spin-precession and cyclotron 
frequencies for polarized muons in a magnetic storage ring. The storage ring magnetic field is measured 
using nuclear magnetic resonance probes calibrated in terms of the equivalent proton spin precession 
frequency @/, in a spherical water sample at 34.7 °C. The ratio w,/@/,, together with known fundamental 
constants, determines a,(FNAL) = 116592040(54) x 1071! (0.46 ppm). The result is 3.3 standard 
deviations greater than the standard model prediction and is in excellent agreement with the previous 
Brookhaven National Laboratory (BNL) E821 measurement. After combination with previous measure- 
ments of both w* and 4^, the new experimental average of a,(Exp) = 116592061(41) x 107" 
(0.35 ppm) increases the tension between experiment and theory to 4.2 standard deviations. 


DOI: 10.1103/PhysRevLett.126.141801 


I. INTRODUCTION 


The magnetic moments of the electron and muon 
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FIG. 1. Feynman diagrams of representative SM contributions 
to the muon anomaly. From left to right: first-order QED and 
weak processes, leading-order hadronic (H) vacuum polarization, 
and hadronic light-by-light contributions. 


the Dirac equation [1] was its prediction for the electron 
that ge = 2. Motivated in part by anomalies in the hyperfine 
structure of hydrogen [2,3], Schwinger [4] proposed an 
additional contribution to the electron magnetic moment 
from a radiative correction, predicting the anomaly [5] 
a, = a[2z c 0.001 16 in agreement with experiment [6]. 

The first muon spin rotation experiment that observed 
parity violation in muon decay [7] determined that, to 
within 1096, g, — 2, which was subsequently measured 
with higher precision [8]. A more precise experiment [9] 
confirmed Schwinger's prediction for the muon anomaly 
and thereby established for the first time the notion that a 
muon behaved like a heavy electron in a magnetic field. 
This evidence, combined with the discovery of the muon 
neutrino [10], pointed to the generational structure of 
the SM. 

The SM contributions to the muon anomaly, as illustrated 
in Fig. 1, include electromagnetic, strong, and weak inter- 
actions that arise from virtual effects involving photons, 
leptons, hadrons, and the W, Z, and Higgs bosons [11]. 
Recently, the international theory community published 
a comprehensive [12-24] SM prediction [13] for the 
muon anomaly, finding a,(SM)—116591810(43) x 107"! 
(0.37 ppm). It is based on state-of-the-art evaluations of the 
contributions from quantum electrodynamics (QED) to 
tenth order [25,26], hadronic vacuum polarization [27— 
33], hadronic light-by-light [11,3447], and electroweak 
processes [48-52]. 

The measurement of a, has become increasingly precise 
through a series of innovations employed by three exper- 
imental campaigns at CERN [53-55] and more recently at 
Brookhaven (BNL E821) [56]. The BNL net result, with its 
0.54 ppm precision, is larger than a,,(SM) by 3.7 standard 
deviations (c). While the electron magnetic anomaly has 
been measured to fractions of a part per billion [57], the 
relative contribution of virtual heavy particles in many 
cases scales as (m, /m,)? ~ 43 000. This is the case e.g. for 
the W and Z bosons of the SM and many hypothetical new 
particles, and it gives the muon anomaly a significant 
advantage when searching for effects of new heavy physics. 
Because the BNL result hints at physics not included in the 
SM, Experiment E989 [58] at Fermilab was constructed to 
independently confirm or refute that finding. In this paper, 
we report our first result with a precision of 0.46 ppm. 


Separate papers provide analysis details on the muon 
precession [59], the beam dynamics corrections [60], 
and the magnetic field [61] determination. 


II. EXPERIMENTAL METHOD 


The experiment follows the BNL concept [56] and uses 
the same 1.45 T superconducting storage ring (SR) magnet 
[62], but it benefits from substantial improvements. These 
include a 2.5 times improved magnetic field intrinsic 
uniformity, detailed beam storage simulations, and state- 
of-the-art tracking, calorimetry, and field metrology for the 
measurement of the beam properties, precession frequency, 
and magnetic field [58]. 

The Fermilab Muon Campus delivers 16 highly polar- 
ized, 3.1 GeV/c, ~120 ns long positive muon beam 
bunches every 1.4 s into the SR. A fast pulsed-kicker 
magnet deflects the muon bunch into a 9-cm-diameter 
storage aperture, resulting in ©5000 stored muons per fill. 
The central orbit has a radius of Rọ = 7.112 m and the 
cyclotron period is 149.2 ns. Four sections of electrostatic 
quadrupole (ESQ) plates provide weak focusing for vertical 
confinement. 

The muon spins precess in the magnetic field at a rate 
greater than the cyclotron frequency. The anomalous 
precession frequency [63] in the presence of the electric 


E and magnetic B fields of the SR is 


8, = 8,- ö, = -L jab- a, (LL) 8-8 


For horizontally circulating muons in a vertical magnetic 


field, B - B = 0; this condition is approximately met in our 
SR. At the muon central momentum pg, set such that 
y, = (1 + 1/a,) © 29.3, the third term vanishes. 

In-vacuum straw tracker stations located at azimuthal 
angle @ = 180° and 270° with respect to the injection point 
provide nondestructive, time-in-fill dependent beam pro- 
files M(x,y,@,t) by extrapolation of decay positron 
trajectories to their upstream radial tangency points within 
the storage aperture [64]. These profiles determine the 
betatron oscillation parameters necessary for beam dynam- 
ics corrections and the precession data fits discussed below. 

Twenty-four calorimeters [65-67], each containing 
a 9 x 6 array of PbF, crystals, detect the inward-spiraling 
decay positrons. When a signal in a silicon photomultiplier 
(SiPM) viewing any crystal exceeds ~50 MeV, the data- 
acquisition system stores the 54 waveforms from that 
calorimeter in a set time window around the event. 
Decay positron hit times and energies are derived from 
reconstruction of the waveforms. 
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The magnetic field is measured using pulsed proton 
nuclear magnetic resonance (NMR), calibrated in terms of 
the equivalent precession frequency w,'(7,) of a proton 
shielded in a spherical sample of water at a reference 
temperature T, = 34.7?C. The magnetic field B is deter- 
mined from the precession frequency and shielded proton 
magnetic moment, jJ, (T,) using io, = 2u,B. The muon 
anomaly can then be obtained from [68-72] 


"E Wa EU) He(H) My Je (2) 
"TEQTAG) a mi 


where our collaboration measures the two quantities to 
form the ratio 


anro) M 


The Run-1 data, collected in 2018, are grouped into four 
subsets (a-d) that are distinguished by unique kicker and 
ESQ voltage combinations. The ratio Rj, can be concep- 
tually written in terms of measured quantities and correc- 
tions as 


R! fase (1 Ce C, Cui Cas) 


x : 4 
M PUCK CET CSC) er ea 59 


The numerator includes the master clock unblinding factor 
Flock» the measured precession frequency œ, and four 
beam-dynamics corrections, C;. We deconstruct @,(T) into 
the absolute NMR calibration procedure (indicated by 
fcaib) and the field maps, which are weighted by the 
detected positrons and the muon distribution averaged over 
several timescales [(@,(x, y, p) x M(x, y, ))]. The result 
must be corrected for two fast magnetic transients B; that 
are synchronized to the injection. 

Damage to two of the 32 ESQ high-voltage resistors was 
discovered after completion of Run-1. This led to slower- 
than-designed charging of one of the quadrupole sections, 
spoiling the symmetry of the electric field early in each fill. 
The impact of this is accounted for in the analysis 
presented. Brief summaries of the terms in Eq. (4) follow. 


III. ANOMALOUS PRECESSION FREQUENCY 


f clocks A single 10 MHz, GPS-disciplined master clock 
drives both the œ, and @/, measurements. The clock has a 
one-week Allan deviation [73] of 1 ppt. Two frequencies 
derived from this clock provide the 61.74 MHz field 
reference and a blinded “(40-e) MHz" used for the c, 
precession measurement. A blinding factor in the range 
+25 ppm was set and monitored by individuals external to 
our collaboration. f sock 1s the unblinding conversion factor; 
its uncertainty is negligible. 


07 —The signature of muon spin precession stems from 
parity violation in 4^ decay, which correlates the muon 
spin and the positron emission directions in the 4^ rest 
frame. When boosted to the lab frame, this correlation 
modulates the e* energy (E) spectrum at the relative 
precession frequency @, between the muon spin and 
momentum directions. The rate of detected positrons 
with E > Ey, as a function of time t into the muon fill 
then varies as 


N(t) = Nony (i)e 
x {1 + Ana (t) cos|ost + po + ng(t)]}, (5) 


where yz, is the time-dilated muon lifetime (~64.4 us), No 
is the normalization, A is the average weak-decay asym- 
metry, and @ is the ensemble average phase angle at 
injection. The latter three parameters all depend on Ep. The 
n; terms model effects from betatron oscillations of the 
beam, and are not required in their absence. This beam 
motion couples with detector acceptance to modulate the 
rate and the average energy, and hence the average 
asymmetry and phase, at specific frequencies. The coherent 
betatron oscillation (CBO) in the radial direction dominates 
the modulation. 

The CBO, aliased vertical width (VW), and vertical 
mean ((y)) frequencies are well measured, and the y; terms 
are well modeled and minimally correlated in fits for @,. 

An accurate fit to the data also requires accounting for 
the continuous loss of muons over a fill, also weakly 
coupled to w,. Coincident minimum-ionizing energies in 
three sequential calorimeters provide a signal to determine 
the time dependence of muon losses. 

Two complementary reconstruction algorithms trans- 
form the digitized SiPM waveforms into positron energies 
and arrival times. In the “local” approach, waveforms are 
template-fit to identify all pulses in each crystal, which are 
then clustered based on a time window. In the “global” 
approach, waveforms in a 3 x 3 array of crystals centered 
on a local maximum in time and position are template-fit 
simultaneously. After subtraction of the fit from the wave- 
forms, that algorithm iterates to test for any missed pulses 
from multiparticle pileup. To avoid biasing w,, we stabilize 
the calorimeter energy measurement within a muon fill by 
correcting the energy reconstruction algorithm on the SiPM 
pixel recovery timescale (up to tens of nanoseconds) and 
the fill timescale (700 us) using a laser-based monitoring 
system [74]. The system also provides long-term (many- 
days) gain corrections. The two reconstructed positron 
samples are used in four independent extractions of œ, in 
which each e* contribution to the time series is weighted 
by its energy-dependent asymmetry; this is the optimal 
approach [75]. Seven other determinations using additional 
methods agree well [59]. Each time series is modified to 
statistically correct for contributions of unresolved pileup 
clusters that result from multiple positrons proximate in 
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FIG. 2. Fourier transform of the residuals from a time-series fit 
following Eq. (5) but neglecting betatron motion and muon loss 
(red dashed), and from the full fit (black). The peaks correspond 
to the neglected betatron frequencies and muon loss. Inset: 
asymmetry-weighted e* time spectrum (black) from the Run- 
1c run group fit with the full fit function (red) overlaid. 


space and time. The analyses employ one of three 
data-driven techniques to correct for pileup, which would 
otherwise bias w,. 

A y? minimization of the data model of Eq. (5) to the 
reconstructed time series determines the measured (m) 
quantity œ. The model fits the data well (see inset to 
Fig. 2), producing reduced ys consistent with unity. 
Fourier transforms of the fit residuals show no unmodeled 
frequency components, see Fig. 2. Without the 7; terms and 
the muon loss function in the model, strong signals emerge 
in the residuals at expected frequencies. 

The dominant systematic uncertainties on @, arise from 
uncertainties in the pileup and gain correction factors, the 
modeling of the functional form of the CBO decoherence, 
and in the @cgo(t) model. Scans varying the fit start and 
stop times and across individual calorimeter stations 
showed no significant variation in any of the four run 
groups [59]. 

The measured frequency œ% requires four corrections, 
C;, for interpretation as the anomalous precession fre- 
quency o, of Eq. (2). The details are found in Ref. [60]. 

C,.—The electric-field correction C, from the last term 
in Eq. (1) depends on the distribution of equilibrium radii 
X, = x — Ro, which translates to the muon beam momen- 
tum distribution via Ap/ po = x,(1— n)/Ro, where n is the 
field index determined by the ESQ voltage [60]. A Fourier 
analysis [60,76] of the decoherence rate of the incoming 
bunched beam as measured by the calorimeters provides 
the momentum distribution and determines the mean 
equilibrium radius (x,) «6 mm and the width c, %9 mm. 
The final correction factor is C, = 2n(1 — n)ff (x2)/ R7, 
where (x2) = o2. + (xey. 

C,,.—A pitch correction C, is required to account for 
the vertical betatron oscillations that lead to a nonzero 


average value of the B - B term in Eq. (1). The expression 


C, = n(A2)/ARj determines the pitch correction factor 
[60,77]. The acceptance-corrected vertical amplitude A, 
distribution in the above expression is measured by the 
trackers. 

Extensive simulations determined the uncertainties ôC, 
and óC, arising from the geometry and alignment of the 
plates, as well as their voltage uncertainties and non- 
linearities. The nonuniform Kicker time profile applied to 
the finite-length incoming muon bunch results in a corre- 
lation introducing the largest uncertainty on C,. 

C,,.—Any bias in the average phase of muons that 
are lost compared to those that remain stored creates 
a time dependence to the phase factor gy in Eq. (5). 
Beamline simulations predict a phase-momentum correla- 
tion dyo/dp — (-10.0+ 1.6) mrad/ (96^ p/po) and losses 
are known to be momentum dependent. We verified the 
correlation by fitting precession data from short runs in 
which the storage ring magnetic field, and thus the central 
stored momentum po, varied by +0.67% compared to its 
nominal setting. Next, we measured the relative rates of 
muon loss (ml) versus momentum in dedicated runs in 
which muon distributions were heavily biased toward high 
or low momenta using upstream collimators. Coupling the 
measured rate of muon loss in Run-1 to these two 
correlation factors determines the correction factor C,,;. 

Cpa The phase term @ in Eq. (5) depends on the 
muon decay coordinate (x, y, p) and positron energy, but 
the precession frequency o, does not. If the stored muon 
average transverse distribution and the detector gains are 
stable throughout a fill, that average phase remains con- 
stant. The two damaged resistors in the ESQ system caused 
slow changes to the muon distribution during the first 
~100 us of the measuring period. An extensive study of 
this effect involved (a) generation of phase, asymmetry, and 
acceptance maps for each calorimeter as a function of muon 
decay coordinate and positron energy from simulations 
utilizing our GEANT-based model of the ring (GM2RINGSIM); 
(b) extraction of the time dependence of the optical lattice 
around the ring from the cosv simulation package and 
GM2RINGSIM; (c) folding the azimuthal beam distribution 
derived from tracker and optics simulations with the phase, 
asymmetry, and acceptance maps to determine a net 
effective phase shift versus time-in-fill, g(t); and (d) appli- 
cation of this time-dependent phase shift to precession data 
fits to determine the phase-acceptance (pa) correction Cpa. 
The use of multiple approaches confirmed the conclusions; 
for details, see Ref. [60]. The damaged resistors were 
replaced after Run-1, which significantly reduces the 
dominant contribution to C,, and the overall magnitude 
of muon losses. 


IV. MAGNETIC FIELD DETERMINATION 


A suite of pulsed-proton NMR probes, each optimized 
for a different function in the analysis chain, measures the 
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magnetic field strength [61]. Every ~3 days during data 
taking, a 17-probe NMR trolley [78] measures the field at 
about 9000 locations in azimuth to provide a set of 2D field 
maps. 378 pulsed-NMR probes, located 7.7 cm above and 
below the storage volume, continuously monitor the field at 
72 azimuthal positions, called stations. The trolley and 
fixed probes use petroleum jelly as an NMR sample. The 
probe signals are digitized and analyzed [79] to extract a 
precession frequency proportional to the average magnetic 
field over the NMR sample volume. A subset of probes is 
used to provide feedback to the magnet power supply to 
stabilize the field. 

Calibration procedure fcap. —The primary calibration 
uses a probe with a cylindrical water sample. Corrections 
are required to relate its frequencies to the precession 
frequency expected from a proton in water at the reference 
temperature 34.7 °C. Studies of the calibration probe in an 
MRI solenoid precisely determine corrections for sample 
shape, temperature, and magnetization of probe materials to 
an uncertainty of 15 ppb. Cross-calibrations to an absolute 
*He magnetometer [80] confirm the corrections to better 
than 38 ppb. 

The calibration probe is installed on a translation stage in 
the SR vacuum. We repeatedly swap the calibration probe 
and a trolley probe into the same location, compensating for 
changes of the SR field. This procedure determines 
calibration offsets between individual trolley probes and 
the equivalent o, values. The offsets are due primarily to 
differences in diamagnetic shielding of protons in water 
versus petroleum jelly, sample shape, and magnetic per- 
turbations from magnetization of the materials used in 
the probes and trolley body. The trolley probe calibration 
offsets are determined with an average uncertainty of 
30 ppb. 

Field tracking [o (x. y, 9) ].—The 14 Run-1 trolley field 
measurements bracket muon storage intervals f, to fy, . 
They provide a suite of 2D multipole moments (dipole, 
normal quadrupole, skew quadrupole, etc.), which the fixed 
probes track. The fixed probes provide five independent 
moments (four moments for some stations) that track the 
field over 5? in azimuth for each station. The trolley 
moments are interpolated for times between the trolley 
runs, and the fixed probes continuously track changes to 
five lower-order moments [61]. The fixed probe and trolley 
measurements are synchronized when the trolley passes, 
averaged over each 5? azimuthal segment. The trolley run at 
time /,,, yields a second set of moments mi (t,.1). 
The fixed probe moments m?(t, p) are used to interpolate 
the field during muon storage between the trolley runs. The 
uncertainty on the interpolation is estimated from both the k 
and k+1 maps and a Brownian bridge random walk 
model. The procedure produces interpolated storage vol- 
ume field maps c(x, y, 9) in terms of the equivalent 
shielded proton frequency throughout the Run-1 data- 
taking periods. 
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FIG.3. Azimuthally averaged magnetic field contours o, (x, y) 
overlaid on the time and azimuthally averaged muon distribution 
M(x, y). 


Muon weighting [M(x, y, $)].—Averaging of the mag- 
netic field weighted by the muon distribution in time and 
space uses the detected positron rates and the muon beam 
distribution measured by the trackers. The interpolated field 
maps are averaged over periods of roughly 10 s and 
weighted by the number of detected positrons during the 
same period. The SR guide fields introduce azimuthal 
dependencies of the muon distribution M(x, y. 9). We 
determine the muon-weighted average magnetic field by 
summing the field moments m; multiplied by the beam- 
weighted projections k; for every three-hour interval over 
which the tracker maps and field maps are averaged. Along 
y, the beam is highly symmetric and centered, and the skew 
field moments (derivatives with respect to y) are relatively 
small. The azimuthally averaged centroid of the beam is 
displaced radially, leading to relative weights for the 
field dipole, normal quadrupole, and normal sextupole of 
k; — 1.0, 0.15, and 0.09, respectively. An overlay of the 
azimuthally averaged field contours on the muon distribu- 
tion is shown in Fig. 3. The combined total uncertainty of 
a’, from probe calibrations, field maps, tracker alignment 
and acceptance, calorimeter acceptance, and beam dynam- 
ics modeling is 56 ppb. 


TABLEI. Run-1 group measurements of @,, à, and their ratios 


Rj, multiplied by 1000. See also Supplemental Material [81]. 


Run @a/27 [Hz] à, /2x [Hz] Ri, x 1000 

la 229 081.06(28) 61791 871.2(7.1) 3.707 300 9(45) 
1b 229 081.40(24) 61791 937.8(7.9) 3.707 302 4(38) 
Ic 229 081.26(19) 61791 845.4(7.7) 3.707 305 7(31) 
ld 229 081.23(16) 61 792003.4(6.6) 3.707 295 7(26) 
Run-1 3.707 300 3(17) 
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TABLE II. Values and uncertainties of the Ri correction terms 
in Eq. (4), and uncertainties due to the constants in Eq. (2) for a,. 


Positive C; increase a, and positive B; decrease a,. 


Correction Uncertainty 
Quantity terms (ppb) (ppb) 
07. (statistical) ee 434 
œ% (systematic) ee 56 
C, 489 53 
C, 180 13 
Cmi -11 5 
Cpa —158 75 
Fei (mp (x. y, d) x M(x, y. d)) UC 56 
B, —27 37 
B, -17 92 
i (34.7°)/pe T 10 
m,/m, ZE 22 
ge/2 e: 0 
Total systematic zn 157 
Total fundamental factors b 25 
Totals 544 462 


B, and B,.—Two fast transients induced by the dynam- 
ics of charging the ESQ system and firing the SR kicker 
magnet slightly influence the actual average field seen by 
the beam compared to its NMR-measured value as 
described above and in Ref. [61]. An eddy current induced 
locally in the vacuum chamber structures by the kicker 
system produces a transient magnetic field in the storage 
volume. A Faraday magnetometer installed between the 
kicker plates measured the rotation of polarized light in a 
terbium-gallium-garnet crystal from the transient field to 
determine the correction B,. 


BNL g-2 —+——_@-——____+- 
FNAL g-2 +————@——+ 


4.20 


——. So] 
Standard Model Experiment 
Average 


175 180 185 190 195 200 205 210 215 
9 
a, * 10 — 1165900 


FIG. 4. From top to bottom: experimental values of a, from 
BNL E821, this measurement, and the combined average. The 
inner tick marks indicate the statistical contribution to the total 
uncertainties. The Muon g — 2 Theory Initiative recommended 
value [13] for the standard model is also shown. 


The second transient arises from charging the ESQs, 
where the Lorentz forces induce mechanical vibrations in 
the plates that generate magnetic perturbations. The ampli- 
tudes and sign of the perturbations vary over the two 
sequences of eight distinct fills that occur in each 1.4 s 
accelerator supercycle. Customized NMR probes measured 
these transient fields at several positions within one ESQ 
and at the center of each of the other ESQs to determine 
the average field throughout the quadrupole volumes. 
Weighting the temporal behavior of the transient fields 
by the muon decay rate, and correcting for the azimuthal 
fractions of the ring coverage, 8.5% and 43% respectively, 
each transient provides final corrections B, and B, to a, as 
listed in Table II. 


V. COMPUTING a, AND CONCLUSIONS 


Table I lists the individual measurements of œ, and &,, 
inclusive of all correction terms in Eq. (4), for the four run 
groups, as well as their ratios, 7v, (the latter multiplied by 
1000). The measurements are largely uncorrelated because 
the run-group uncertainties are dominated by the statistical 
uncertainty on o. However, most systematic uncertainties 
for both œ, and à", measurements, and hence for the ratios 
R are fully correlated across run groups. The net computed 
uncertainties (and corrections) are listed in Table II. The fit 
of the four run-group results has a y*/n.d.f. = 6.8/3, 
corresponding to P(y”) = 7.8%; we consider the P(x?) to 
be a plausible statistical outcome and not indicative of 
incorrectly estimated uncertainties. The weighted-average 
value is R, = 0.003 707 3003(16) (6), where the first error 
is statistical and the second is systematic [82]. From Eq. (2), 
we arrive at a determination of the muon anomaly 

a,(FNAL) = 116592040(54) x 10"!! (0.46 ppm), 
where the statistical, systematic, and fundamental constant 
uncertainties that are listed in Table II are combined in 
quadrature. Our result differs from the SM value by 3.36 and 
agrees with the BNL E821 result. The combined exper- 
imental (Exp) average [83] is 


a,(Exp) = 116592061(41) x 107!! (0.35 ppm). 


The difference, a, (Exp) — a,(SM) = (251 +59) x 107!!, 
has a significance of 4.20. These results are displayed 
in Fig. 4. 

In summary, the findings here confirm the BNL exper- 
imental result and the corresponding experimental average 
increases the significance of the discrepancy between the 
measured and SM predicted a, to 4.20. This result will 
further motivate the development of SM extensions, 
including those having new couplings to leptons. 

Following the Run-1 measurements, improvements to 
the temperature in the experimental hall have led to greater 
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magnetic field and detector gain stability. An upgrade to 
the kicker enables the incoming beam to be stored in the 
center of the storage aperture, thus reducing various beam 
dynamics effects. These changes, amongst others, will lead 
to higher precision in future publications. 
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